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I.  INTRODUCTION 


1 


A.  GENERAL  FEATURES  OF  THE  ANOMALOUS  PHOTOVOLTAIC  EFFECT 

Certain  crystals  of  ZnS  generate  voltages  greater  than  the  band  gap  voltage 
when  illuminated  with  ultraviolet  light.  This  phenomenon  is  called  the  anoma¬ 
lous  photovoltaic  effect.  ^ 

2  3 

It  is  well  known,  from  thermodynamic  arguments,  ’  that  the  maximum  photo¬ 
voltage  a  semiconductor  device  can  generate  is  equal  to  its  band  gap  voltage. 

4 

Further  studies  by  Merz  show  that  this  photovoltage  can  be  "quenched”  by 
simultaneously  illuminating  the  crystal  with  more  than  one  wavelength  of  light. 
Cheroff  and  Keller  found  that  the  polarity  of  the  photovoltage  is  dependent  on 
the  wavelength  of  the  incident  light.  Figure  1  is  a  graph  of  the  open-circuit 
voltage  and  short-circuit  current  as  a  function  of  the  wavelength  of  the  incident 
light,  for  a  typical  crystal  used  in  the  present  study. 

In  general,  it  is  found  that  at  wavelengths  longer  than  350  m/x,  the  photovoltages 
are  of  the  order  of  1  volt  or  less .  At  about  350  mp  the  photovoltage  reverses 
polarity  and  reaches  a  maximum  at  approximately  340  xafi  .  This  maximum  (to 
be  called  the  positive  peak)  can  be  of  the  order  of  100  volts.  At  shorter  wave¬ 
lengths  (approximately  330  m/x)  the  photovoltage  again  changes  polarity  and  in¬ 
creases  to  a  maximum  (of  the  opposite  polarity).  This  negative  peak  is  often 
of  the  same  order  of  magnitude  as  the  positive  peak,  but  can  be  greater  in  abso¬ 
lute  magnitude  than  the  positive  peak. 

At  still  shorter  wavelengths,  the  behavior  varies.  In  some  crystals,  the  photo¬ 
voltage  decreases  to  zero  and  in  others  the  photovoltage  persists  to  very  short 
wavelengths  (approximately  280  m/x ) . 
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Fig.  1  Spectral  Dependence  of  Anomalous  Photovoltaic  Effect 


Typical  spectral  dependence  of  the  open-circuit  voltage  Vqc  ^nd  short-circuit 
current  for  a  ZnS  crystal  exhibiting  anomalous  photovoltaic  effect. 
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The  purpose  of  this  study  will  be  to  devise  and  justify  a  naodel  (semi-quantitative) 
to  account  for  the  general  features  of  the  anomalous  photovoltaic  effect. 

B.  SOME  PROPERTIES  OF  ZnS 

Before  reviewing  more  of  the  work  done  on  the  anomalous  photovoltaic  effect 
(APE),  a  survey  of  some  of  the  properties  of  zinc  sulfide  will  be  helpful. 

Zinc  sidfide  exists  in  both  the  cubic  (sphalerite)  and  hexagonal  (wurtzite)  struc- 
tures.  Allen  and  Crenshaw  find  that  the  cubic  phase  transforms  into  the  hexa¬ 
gonal  phase  on  heating  to  1020 "C.  This  figure  for  the  transition  temperature 

must  be  viewed  with  some  reservation,  since  it  has  been  found  that  the  transition 

7 

temperature  is  dependent  on  the  purity  of  the  sample  used.  Allen  and  Crenshaw's 
work  was  done  on  a  sphalerite  specimen  containing  0. 15%  Fe.  The  presence  of 
chlorides  either  lowers  the  transition  temperature  or  at  least  speeds  the  transi- 

7 

tion.  Kremheller  finds  no  transformation  even  at  temperatures  as  high  as 

1175°C  in  the  absence  of  chlorine,  a  common  contaminant  in  even  the  purest  ma- 

8  9 

terial.  There  is  also  some  evidence  ’  that  the  hexagonal  phase  transforms 
back  to  the  cubic  phase  at  1240°C.  One  important  feature  of  the  transformation 
at  1020 °C  is  that  the  transformation  back  to  cubic  from  hexagonal  on  cooling 
from  above  1020°C  is  extremely  slow.  In  fact,  it  is  possible  to  cool  the  sample 
to  room  temperature  without  transformation  into  the  cubic  phase. 

It  is  significant  that  only  crystals  grown  from  the  vapor  phase  show  anomalous 
behavior.  This  method  generally  produces  needle-like  crystals  up  to  1  cm  long. 
The  diameter  of  the  crystals  is  usually  a  fraction  of  a  millimeter.  Structurally 
perfect  crystals  (either  hexagonal  or  cubic)  are  rarely,  if  ever,  produced  by  this 

g 

method.  Usually  the  crystals  display  some  kind  of  stacking  disorder.  The  one 
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characteristic  all  crystals  grown  in  this  manner  have  in  common  is  that  the  long 
axis  of  the  crystal  is  either  the  cubic  [ill]  axis  or  the  hexagonal  "c"  axis. 

Some  information  about  the  optical  band  edge  of  both  phases  of  ZnS  is  essential 
to  the  analysis  of  the  APE.  Krdger^^  finds  that  cubic  material  displays  a  sharp 
edge  at  340  m//  .  This  result  is  confirmed  by  Beun  and  Goldsmith.  The  re¬ 
sults  of  measurements  of  the  optical  band  edge  of  hexagonal  material  by  vari¬ 
ous  researchers  are  not  in  as  good  agreement  as  the  measurements  on  cubic 

5  12  13  14 

material.  Piper  and  others  *  ’  ’  find  an  edge  at  approximately  335  mfi . 

However,  Beun  and  Goldsmith^^  place  the  edge  at  approximately  330  mfi.  Since 
structurally  pure  crystals  are  rarely  found,  and  the  presence  of  any  cubic  ma¬ 
terial  will  shift  the  edge  to  longer  wavelengths,  it  is  felt  that  the  results  of 

Beun  and  Goldsmith  are  the  most  reliable .  The  effect  of  polarized  light  on  the 

11  14 

position  of  the  band  edge  is  also  known  and  will  be  useful.  ’ 

C.  ESSENTIAL  FEATURES  OF  THE  APE 

4 

The  most  significant  observation  to  come  out  of  the  early  research  is  the  fact 

that  only  crystals  displaying  "banding"  (when  viewed  between  crossed  polarizers) 

produce  anomalous  effects.  Such  crystals  are  thought  to  be  composed  of  layers 

15  16  IT 

of  hexagonal  and  cubic  material.  ’  A.  Lempicki  finds,  by  means  of  x  ray 

analysis,  that  pure  hexagonal  crystals  show  no  anomalies.  On  the  evidence  siq)- 

4 

plied  by  his  own  optical  observations,  Merz  proposes  that  the  anomalously  high 
photovoltages  can  be  accounted  for  on  a  model  which  postulates  the  additivity  of 
voltages  generated  by  some  simple,  normally  behaving  unit.  Thus,  by  counting 

4 

the  bands  observed  on  crystals  between  crossed  polarizers,  Merz  estimates 
that  the  voltage  due  to  one  unit  must  be  approximately  0. 1  volt.  Furthermore, 
if  a  crystal  which  displays  uniform  banding  is  broken  in  half,  each  part  produces 
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approximately  half  the  photovoltage  observed  for  the  whole  crystal.*  Merz's 

4 

own  model,  which  does  not  account  for  the  spectral  response  curve  (it  had  not 

been  measured  at  that  time),  cannot  account  for  the  high  voltages  since,  as 
18 

pointed  out  by  Tauc ,  it  predicts  zero  photovoltage  at  high  light  levels .  This  is 

in  contradiction  to  Merz's  own  data. 

4 

Another  significant  observation  is  that  high  voltages  are  recorded  only  along  the 

long  ("c")  axis  of  the  crystal.  This  is  confirmed  by  measurements  on  platelets 

of  ZnS.  For  these  crystals,  the  "c"  axis  can  be  determined  optically.  Merz 

also  noted  that  there  is  a  definite  relation  between  the  direction  of  growth  and 

the  polarity  of  the  photovoltage  (in  his  work,  this  is  the  positive  peak  photovolt- 
17 

age) .  liCmpicki  finds  no  such  correlation,  but  notes  that  the  direction  of  the 
positive  peak  is  correlated  with  the  direction  of  the  "c"  axis  as  determined  by 
pyroelectric  raeasuxements.  In  this  ejqxeriment,  Lempicki  found  that  the  end  of 
the  crystal  which  became  positive  on  heating  was  positive  at  340  tnfi .  All  re¬ 
searchers^’  find  that  the  positive  peak  falls  at  a  wavelength  near  the 

absorption  edge  previously  reported  for  hexagonal  ZnS. 


17 

♦Lempicki  ,  however,  states  that  merely  picking  up  a  crystal  with  forceps 
can  alter  its  characteristics.  Thus,  Merz's  results  on  breaking  crystals  could 
be  questioned. 
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II.  EXPERIMENTAL  PROCEDURE 

A.  SAMPLE  PREPARATION 

ZnS  crystals  can  be  grown  by  the  sealed  tube  sublimation  technique  described 

13 

by  Reynolds  and  Czyzak.  RCA  "Luminescent"  grade  ZnS  powder,  the  purest 
commercially  available,  is  the  starting  material.  This  powder  is  free  from  acti¬ 
vators.  The  ZnS  powder  is  placed  in  a  clean,  dry,  quartz  tube  of  the  form  shown 

-5 

in  Fig.  2a-  The  tube  is  evacuated  to  a  pressure  of  ^proximately  10  Torr 
and  both  the  tube  and  ZnS  powder  outgassed  by  heating  to  about  900°C .  The  tube 
is  then  back-filled  with  about  10  Torr  of  H  S,  sealed  off  and  placed  in  a  tubular 
furnace  possessing  a  temperature  gradient.  The  charge  of  ZnS  is  placed  in  the 
hottest  region  of  the  furnace  and  the  temperature  there  maintained  at  about 
1170°C.  At  this  temperature,  the  vapor  pressure  of  ZnS  is  nearly  10  Torr 
(Ref.  19).  The  ZnS  vapor  condenses  on  the  coolest  surfaces  of  the  tube,  and,  under 
the  proper  conditions,  crystals  are  formed  at  about  1000°  C.  It  normally  takes  120 
hours  to  produce  crystals  of  usable  dimensions.  The  role  played  by  the  HgS  is  not 
known,  but  it  does  not  aet  as  a  chemical  transport  gas  as  does  iodine  in  the  work 
of  Nilsche,  since  the  material  transport  rate  decreases  with  increasing  HgS 
pressure. 

The  peeuliar  shape  of  the  receiving  surfaces  (Fig.  2a)  enhances  crystal  growth. 
Attempts  to  grow  crystals  in  tubes  without  the  short  inner  tube  (Fig.  2b)  are 
much  less  successful.  Instead  of  obtaining  single  crystals,  a  heavy  layer  of 
polycrystalline  material  is  deposited  on  surface  C.  In  the  type  of  tube  shown  in 
Fig.  2a,  crystals  appear  on  surface  A,  and  a  thin  polycrystalline  layer  forms 
on  surface  B.  Surfaces  B  and  C  are  cooler  than  surface  A,  since  the  former 


LOCKHEED  MISSILES  ft  SPACE  COMPANY 


6-90-63-20 


7 


Fig.  2  Crystal  Growing  Tubes 


Typical  tubes  used  in  the  sublimation  method  of  growing  ZnS  crystals.  Tubes  of 
tj^  "a"  produced  better  crystals. 
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"see"  the  cold  end  of  the  furnace  and  can  be  cooled  by  radiation.  Surface  A 
radiates  radially  toward  the  1000 °C  inner  surface  of  the  furnace.  The  fact  that 
surfaces  B  and  C  are  cool  decreases  the  allowable  vapor  pressure  locally,  thus 
enhancing  the  chance  of  growth.  The  fact  that  surface  A  is  warmer  favors  the 
growth  of  single  crystals. 

B.  MEASUREMENTS 

The  principal  measurement  made  in  this  study  was  the  spectral  response  of  the 
anomalous  photovoltaic  effect.  This  measui'emenl  was  made  at  varying  light 
levels  and  on  both  stressed  and  unstressed  crystals.  Other  auxiliary  measure¬ 
ments  (photoconductive  effect  and  piezoelectric  effect)  will  be  helpful  in  devising 
the  model. 

1 .  Optical  Apparatus 

The  essential  elements  in  the  experimental  equipment  are  a  source  of  mono¬ 
chromatic  light,  aixl  an  electrometer.  A  Leiss  single -pass,  quartz  prism, 
monochromator  in  conjunction  with  a  high-pressure,  water-cooled  mercury  arc 
provides  the  monochromatic  light.  The  lamp  is  rated  at  ^proximately  1000 
watts.  However,  it  is  not  a  continuous  source.  The  lines  are  pressure  broadened 
so  that  there  is  a  continuous  background,  but  there  still  are  some  large  peaks 
in  the  output  of  the  lamp.  As  equipment  is  not  available  to  calibrate  the  lamp 
absolutely,  the  curves  of  spectral  response  depicted  in  the  present  study  are 
uncorrected  for  variations  in  light  intensity. 

The  monochromator  dial  setting  (which  determines  the  output  wavelength)  can  be 
calibrated  by  using  a  standard  spectral  lamp  as  a  line  source,  and  comparing 
the  observed  spectrum  with  a  map  of  the  spectrum  (e.  g. ,  Eder:  Atlas  Typischer 
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Spektren) .  Variations  in  intensity  of  the  light  incident  on  the  sample  are  accom¬ 
plished  by  placing  copper  screens  in  the  light  beam. 

12 

As  the  impedance  of  these  ZnS  crystals  is  often  of  the  order  of  10  ohms,  it  is 

necessary  to  use  an  electrometer  to  measure  the  photovoltage.  The  electrom- 

14 

eter  (Keithly  200B)  has  an  input  impedance  of  approximately  5  x  10  ohms. 

The  short-circuit  current  is  observed  by  measuring  the  voltage  drop  across  a 
calibrated  shunt  resistor.  The  shunt  resistance  must  be  considerably  lower 
than  the  sample  resistance. 

2.  Electrical  Contacts  and  Connections 

Electrical  contact  to  ZnS  can  be  made  with  silver  paint.  These  contacts  give 
the  same  results  as  Hg-In-Tf  liquid  contacts,  and  In-Hg  alloy  solder  contacts. 
The  silver  paint  contacts  are  the  easiest  to  apply  and  were  used  on  almost  all 
of  the  samples . 

A  metal  box,  with  accessories  such  as  the  sample  holder  (Fig.  3,  4)  can  be 
connected  directly  to  the  electrometer  by  means  oi  bayonet-type  plugs.  This 
eliminates  all  leads  outside  the  chassis. 

3.  Stress  Measurements 

Measurements  of  the  spectral  response  of  the  photovoltaic  effect  on  stressed 
crystals  and  measurement  of  the  piezoelectric  constant  "g"  are  made  with  the 
same  apparatus.  Again  the  £q)paratus  is  built  into  a  chassis  which  is  connected 
directly  to  the  electrometer.  The  stressing  apparatus  is  shown  in  Fig.  5  and  is 
self-explanatory.  Silver  paint  is  applied  to  the  crystal  and  the  crystal  placed, 
end  on,  on  the  stainless  steel  pedestal.  The  set  screw  is  then  adjusted  until 
contact  is  made  with  terminal  A.  A  known  mass  is  then  lowered  onto  the 
phosphor  bronze  strip  A'  and  the  desired  measurement  is  made. 
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Fig.  3  Crystal  Support 


Crystals  were  mounted  in  this  manner  for  measurements  of  photovoltaic  effects, 
current-voltage  characteristics,  and  photoconductance. 
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A  -  TEFLON  BODY 
B-  COPPER  SLAB 
C-  PHOSPHOR  BRONZE 
SPRING  FINGERS 
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Fig.  5  Stress  Apparatus  (Longitudinal  Section) 


A, A'-  PHOSPHOR  BRONZE 
B-  TEFLON 

C  -  iTAI  NLESS  STEEL  PEDESTAL 
0  -  BRASS  BASE 
M-  WEIGHT 
S  -  SPECIMEN 


Schematic  drawing  of  apparatus  used  to  apply  stress  to  crystals  in  measurements 
of  the  stress  dependence  of  the  anomalous  photovoltaic  effect  and  the  peizoelectric 
constant  "g." 
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4 .  Measurement  of  Cur  rent- Voltage  Characteristics  and  Photoconductance 

Both  the  current-voltage  and  photoconductance  data  are  taken  with  the  same 
circuit  (Fig.  6).  It  should  be  noted  that  Rg  must  be  kept  approximately  two 
orders  of  magnitude  lower  than  the  sample  resistance.  Rg  is  always  at  least 
an  order  of  magnitude  lower  than  Rg. 
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Fig.  6  Measuring  Circuit 


E  -  ELECTROMETER 

S  -  SPECIMEN 

Rg-  SHUNT  RESISTOR 

Rg-  VOLTAGE  DIVIDER  RESISTANCE 

t  -  APPLIED  VOLTAGE 


Circuit  used  to  measure  current-voltage  characteristics  and  photoconductance 
of  ZnS  crystals. 


LOCKHEED  MISSILES  &  SPACE  COMPANY 


I 

1  m.  PROPOSED  ORIGINS  OF  ANOMALOUS  PHOTOVOLTAIC  EFFECT 


A.  REVIEW  AND  CRITICISM  OF  PROPOSED  MODELS 

Two  postulates  are  common  to  the  three  models  proposed  to  account  for  the 
anomalous  photovoltaic  effect  In  ZnS:  First,  both  hexagonal  and  cubic  ma¬ 
terial  must  be  present  in  the  crystals  in  alternating  layers  {perpendicular  to  the 
"c"  axis);  second,  the  observed  high  voltage  is  the  sum  of  voltages  generated 
in  each  layer  or  between  layers  (called  segments) . 
i8 

Tauc's  model  makes  use  of  the  fact  that  there  is  about  a  0. 1  ev  difference  in 

band  gaps  between  the  hexagonal  and  cubic  phase  of  ZnS  (3.75  ev  and  3.65  ev, 

respectively^^).  Tauc  also  assumes  that  an  impurity  gradient  (donor  or  acceptor) 

21  22 

exists  in  each  segment.  Tauc  shows  ’  that  changes  in  band  gap  and  impurity 
gradients  can  produce  photovoltages.  The  model  is  so  constructed  that  the  im¬ 
purity  gradient  photovoltages  (bulk  photovoltages)  are  opposed  to  the  photovoltage 
generated  by  the  differences  in  band  gap.  The  spectral  response  of  the  APE  can 
be  accounted  for  on  this  model  because  of  the  nature  of  the  band  gap  photovoltage . 
This  voltage  is  generated  only  when  both  the  hexagonal  and  cubic  phases  are 
absorbing. 

On  the  other  hand,  the  bulk  photovoltages  depend  only  on  the  absorption  in  their 

respective  phases.  Thus,  at  long  wavelengths  (e.g. ,  400  m/x )  a  bulk  photovoltage 

in  the  cubic  phase  is  generated.  For  an  impurity  gradient  which  changes  the 

resistivity  of  the  material  by  50%  over  the  length  of  the  segment,  this  voltage  will 
-3  22 

be  about  1.0x10  volts.  As  the  wavelength  is  decreased  to  a  point  where 
the  hexagonal  phase  absorbs,  the  band  gap  photovoltage  is  generated,  and  for 
still  shorter  wavelengths ,  the  hexagonal  phase  becomes  very  absorbing  and  a 
bulk  photovoltage  in  that  phase  is  generated.  Since  the  bulk  and  band  gap  voltages 
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are  opposed,  there  will  be  a  wavelength  at  which  the  band  gap  photovoltage 
(0. 1  ev)  counteracts  the  cubic  bulk  photovoltage  and  the  resultant  voltage  is 
zero.  The  band  gap  photovoltage  then  increases  and  reaches  a  maximum.  The 
maximum  observed  photovoltage  is  the  band  gap  photovoltage  reduced  by  the  two 
bulk  photovoltages . 

In  order  for  the  polarity  of  the  photovoltage  to  change  again,  Tauc  assumed  that 
the  band  gap  photovoltage  decreases  sharply  at  about  330  m^ .  He  proposes  that 
the  electrons  excited  at  this  wavelength  have  sufficient  energy  to  tunnel  through 
the  barrier.  At  very  short  wavelengths  (e.g. ,  <  330  mp)  the  observed  photo¬ 
voltage  is  due  to  the  two  bulk  photovoltages .  The  high  voltage  observed  at  the 
positive  peak  is  then  the  sum  of  the  band  g^  voltage  generated  between  each 
layer. 

This  model  has  one  obvious  drawback.  Essentially,  the  high  positive  peak  photo¬ 
voltage  is  due  to  the  0.1  volt  generated  at  each  band  gap  change  (this  is  the  maximum 

oo  Sect  4  2 

voltage  one  can  generate  from  this  configuration  ’  ‘  *  )•  Since  this  voltage 

is  due  to  the  change  in  band  gap,  it  is  not  reversible  in  the  sense  that  it  also 
could  account  for  the  high  voltages  at  the  negative  peak.  Another  drawback  of 
this  model  is  the  method  of  addition  of  photovoltages .  To  examine  this  facet  of 
the  model,  consider  the  following  band  diagram  for  a  region  of  a  segmented  crystal. 


* 


H  C  H  C  H 


1  I - 1 _ r 


1  1 

A 

B 

■o  X  - ►DISTANCE  ALONG  "c"AXIS 
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The  brdcen  lines  represent  the  impurity  gradients.  At  first  glance  it  appears 
that  the  photovoltage  would  be  exactly  counteracted  by  Vgg,  and  that  no 
addition  could  take  place.  However,  since  the  impurity  density  changes  along  a 
segment,  Aff/ff  (the  relative  change  in  conductivity  under  illumination)  is  also  a 
function  of  x  and  the  combination  of  change  in  band  gap  and  change  in  A  or/o  allows 
one  to  add  the  voltages  (see  Sect.  in-C)  at  each  band  gap  change.  Unfortunately, 
in  order  to  obtain  the  full  band  gap  voltage,  Atr/a  must  vary  considerably  along 
the  segment.  In  fact,  for  a  50%  change  in  conductivity  along  each  segment  (about 
what  is  necessary  to  generate  the  long  wavelength  negative  voltage) ,  a  photo- 
voltage  of  only  (0.  3  AEg)/e  can  be  realized  between  points  A  and  D  -  not  nearly 
enou^  to  support  Merz's  data. 

23 

G.  Neumark  has  proposed  a  model  which  also  uses  the  band  gap  voltage  as  the 
origin  of  the  anomalously  high  positive  peak  photovoltage.  Again,  this  model 
cannot  account  for  the  high  negative  peak  photovoltage.  Addition  of  the  individual 
band  gap  voltages  is  taken  care  of  by  modulating  the  conductivity  by  means  of 
internal  polarization  fields. 

Thus,  both  these  models  break  down  for  the  same  reason.  It  appears  that  any 
model  based  on  the  band  gap  photovoltage  cannot  account  for  all  aspects  of  the 
APE. 

24 

The  final  model  to  be  discussed  was  proposed  by  A.  R.  Hutson.  This  model 
makes  use  of  the  fact  that  both  phases  of  ZnS  are  piezoelectric*  and  that  the  trans¬ 
verse  lattice  constant  of  the  hexagonal  phase  is  smaller  than  the  cubic  lattice 
constant. 


*Cubic  ZnS  is  piezoelectric  only  in  shear,  which  is  all  that  is  needed  in  this  case. 
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Ik  brief,  Hutson  proposes  that  when  there  is  a  crystal  phase  transformation 

during  the  growth  process,  a  transverse  strain  is  frozen  into  the  crystal.  The 

material  on  the  hexagonal  side  of  the  transition  is  in  tension  and  the  material 

on  the  cubic  side  is  in  compression  (see  Fig.  7).  Since  ZnS  is  piezoelectric, 

there  will  be  fields  generated  within  each  segment.  In  the  dark,  these  fields 

are  compensated  for  by  fields  set  up  by  the  diffusion  of  carriers,  as  is  the  case 

22 

in  any  solid  state  photovoltaic  device.  Under  illumination,  however,  these 
fields  can  give  rise  to  a  photovoltage.  Hutson  proposes  that  the  spectral  re¬ 
sponse  be  accounted  for  in  some  way  with  the  difference  in  absorption  coefficients 
between  the  two  phases .  Hutson's  presentation  is  descriptive  and  he  makes  no 
attempt  to  calculate  the  photovoltage  on  the  basis  of  his  model.  He  is  also  not 
specific  on  how  the  photovoltages  add.  In  the  following  sections,  Hutson's 
model  will  be  extended  to  show  how  the  spectral  response  curve  can  be  accounted 
for.  In  addition,  an  expression  for  the  photovoltage  generated  by  a  strained  cry¬ 
stal  will  be  calculated  and  a  model  set  up  to  show  how  the  photovoltages  in  each 
segment  add. 

B.  EXTENSION  OF  HUTSON'S  MODEL 
1 .  Discussion 

The  data  to  be  presented  later  will  be  used  to  justify  a  model  based  on  Hutson's 
ideas.  Before  presenting  the  data,  this  extended  model  will  be  considered  in 
detail . 

This  model  is  predicated  on  two  assumptions: 

(1)  Strains  exist  at  the  transition  between  the  hexagonal  and  cubic  segments. 

(2)  The  hexagonal  segments  are  slightly  more  absorbing  than  the  cubic 
segments  at  wavelengths  of  400  mp  and  longer. 
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Fig.  7  Variation  of  Lattice  Constant  in  Vicinity  of  Structure  Change 


Schematic  variation  of  transverse  lattice  constant  along  "c"  axis  of  a  mixed 
structure  ZnS  crystal,  in  the  vicinity  of  a  structural  change.  This  is  the  type  of 
variation  used  in  the  calculation  of  the  photovoltage  for  the  extended  Hutson 
model. 
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For  simplicity,  the  crystal  will  be  considered  as  composed  of  one  layer  each  of 
cubic  aiKl  hexagonal  material.  Then,  by  Assumption  (1)  and  the  discussion  in  the 
previous  section,  opposing  fields  can  exist  in  the  two  segments .  It  is  these 
fields  which,  when  non-equilibrium  electron-hole  pairs  are  created  by  light, 
give  rise  to  the  photo  voltages.  The  magnitude  of  these  photovoltages  is  depen¬ 
dent  on  the  number  of  electron-hole  pairs  created  (see  Sect.  III-B.  3)  and  thus 
on  the  absorption  coefficients  of  the  material  (Fig.  8).  Since  by  Assumption  (2), 
hexagonal  material  absorbs  more  at  long  wavelengths,  a  small  voltage  will  be 
generated  by  the  field  in  the  hexagonal  segment.  This  voltage  is  chosen  as 
negative  and  is  observed  at  wavelengths  longer  than  350  m/t .  As  the  wavelength 
is  decreased,  the  absorption  edge  of  the  cubic  phase  is  approached,  and  the  cu¬ 
bic  phase  begins  to  absorb.  The  field  in  this  phase  produces  a  photovoltage  which 
counteracts  the  long  wavelength  voltage .  The  resultant  voltage  then  goes  to  zero 
( ~348  m^).  This  wavelength  is  very  close  to  the  region  where  the  absorption 
coefficient  of  cubic  ZnS  increases.  As  the  cubic  material  absorbs  more  heavily, 
its  photovoltage  increases.  This  photovoltage  (denoted  the  positive  peak)  reaches 
its  maximum  at  ~340  m^  (the  cubic  band  edge^^),  then  begins  to  decrease.* 

In  the  region  between  330  and  335  mu,  hexagonal  ZnS  begins  to  absorb.  The  re¬ 
sultant  increase  in  the  voltage  in  this  segment  counteracts  the  positive  peak 
voltage  and  is  the  cause  of  the  decrease  in  observed  photovoltage  in  this  region. 
Eventually  the  voltage  generated  by  the  hexagonal  segment  exactly  counteracts 
the  voltage  in  the  cubic  segment  and  the  observed  photovoltage  is  zero  ( ~330  m/x). 
At  shorter  wavelengths  it  is  maintained  that  the  absorption  coefficient  of  cubic 


♦The  fact  that  the  short-circuit  current  peaks  at  335  mp  is  due  to  a  bright  line 
in  the  Hg  source. 
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Fig.  8  Wavelength  Dependence  of  Transmission  of  Hexagonal  and  Cubic 
ZnS 


Change  in  optical  transmission  in  the  vicinity  of  the  band  edges  of  hexag¬ 
onal  and  cubic  ZnS.  Curve  "a"  is  for  hexagonal  material  with  the  elec¬ 
tric  vector  E  of  the  incident  light  parallel  to  the  "c"  axis,  curve  "b"  is 
for  E±  (c),  and  curve  "d"  is  for  cubic  material  (after  Beun  and  Goldsmith, 
Helv.  Phys.  Acta  3^,  508  (I960)) . 
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ZnS  is  so  large  that  carriers  are  only  created  on  the  surface  and  do  not  participate 
in  the  photovoltaic  effect.  As  the  wavelength  is  decreased,  the  photovoltage  gen¬ 
erated  by  the  hexagonal  segment  reaches  its  maximum  (the  negative  peak).  At 
still  shorter  wavelengths ,  the  absorption  coefficients  of  both  phases  are  so  high 
that  no  carriers  are  generated  in  the  bulk  of  either  phase  and  the  photovoltage 
decreases  to  zero. 

This  argument  can  account  for  the  spectral  response  of  the  APE.  In  order  to 
produce  the  high  voltages  observed  at  the  positive  and  negative  peaks,  note  that  in 
these  two  regions  the  voluige  is  associated  with  one  or  the  other  phase,  not  both. 
That  is,  at  the  positive  peak,  virtually  only  cubic  material  absorbs;  at  the  nega¬ 
tive  peak,  only  the  hexagonal  phase  absorbs  (absorbing  in  the  sense  of  producing 
electron-hole  pairs  in  the  bulk  of  the  crystal).  Thus,  when  conditions  are  right 
for  a  peak,  the  crystal  as  a  whole  is  composed  of  alternating  segments  of  ab¬ 
sorbing  and  non-absorbing  layers.  This  configuration  allows  the  voltages  of 
the  absorbing  segments  to  add  (see  Sect.  m-C). 

One  serious  objection  against  the  model  as  presented  here  is  the  matter  of  sym¬ 
metry  of  strains  in  the  crystal.  That  is,  since  the  crystal  is  composed  of  alter¬ 
nating  layers  of  hexagonal  and  cubic  material,  is  not  the  contribution  to  the 
photovoltage  of  the  strains  at  a  cubic  to  hexagonal  transition  exactly  counteracted 

by  the  contribution  from  a  hexagonal  to  cubic  transition?  Fortimately,  this  is 

0 

not  the  case.  As  pointed  out  by  Allen  and  Crenshaw,  the  hexagonal  to  cubic 
transition  is  very  sluggish,  so  that  the  strain  may  be  spread  over  a  much  larger 
distance  along  the  crystal  at  the  hexagonal  to  cubic  transition.  This  sluggishness 
may  also  allow  impurities  to  migrate  to  this  transition  region  and  release  the 
strain.  Furthermore,  as  will  be  seijn  in  the  next  section,  the  photovoltage  de¬ 
pends  on  the  strain  gradient  rather  than  the  strain.  Thus,  it  is  possible  for 


LOCKHEED  MISSILES  ft  SPACE  COMPANY 


23 


the  cubic  hexagonal  transitions  to  be  sisymmetric  (as  far  as  the  photovoltaic 
effect  is  concerned) . 

2.  Extension  of  Hutson's  Model  —  Derivation  of  Photovoltage 

a.  Method.  The  photovoltage  generated  by  a  nonuniformly  strained 
piezoelectric  crystal  will  be  derived  for  a  one -dimensional  crystal  in  which  the 
strain  falls  off  linearly  with  distance  and  is  zero  for  most  of  the  crystal  (Fig.  9). 
This  photovoltage  will  be  associated  with  one  of  the  segments  of  a  ZnS  crystal. 
Similar  expressions  apply  to  the  other  segments  and  the  observed  photovoltage 
will  be  obtained  by  addition  of  these  elemental  voltages .  The  program  of  the 
derivation  is  as  follows; 

(1)  The  electrostatic  potential  V  inside  the  crystal  will  be  calculated. 

(2)  Since  the  electrochemical  potential  /x  is  a  constant  throughout 
the  segment,  the  slope  of  the  chemical  potential  df /dx  can  be 
found  from  step  (1)  and  the  relations  p  =  f  -eV  and  dfi  =  0. 

(3)  dS/dx  will  be  expressed  in  terms  of  dn/dx,  the  slope  of  the  carrier 
concentration  along  the  crystal. 

21  22 

(4)  Finally,  di./dx  will  be  used  to  calculate  the  photovoltage.  ’ 

b.  Calculation  of  the  Electrostatic  Potential  V.  The  relation  between 
field  E,  displacement  D,  dielectric  constant  e,  piezoelectric  constant  g,  and 
stress  n  is  given  by  Mason;^^ 


E  =  2  _  g  n  (1) 

n  is  related  to  the  strain  S  by  the  relation 

n  =  Ss  , 
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Fig,  9  Variation  of  Strain  Along  Crystal  Segment 


DISTANCE  THROUGH  A  SEGMENT,  i - ► 


Proposed  variation  of  transverse  strain  S  with  distance  along  "c"  axis  of  a  mixed 
crystal  of  ZnS.  This  configuration  was  used  in  the  calculation  of  the  photovoltage 
in  the  extended  Hutson  model. 
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where  s  is  a  constant  (the  elastic  stiffness). 


dE  ^  dS 


Thus,  Poisson's  equation  becomes 
=  4  ir  p  ,  (2) 


where  p  is  the  net  charge  density.  For  an  n-type  semiconductor  with  carrier 

27 

concentration  n  and  N  fully  ionized  donors  per  unit  volume, 

p  =  -  e  (n  -  N)  , 
and  Equation  (2)  can  be  rewritten  as 


dS 


4ire 


(n  -  N) 


Fur  tliermore,  n  is  related  to  £  by 

n  =  N  exp  [^/kX]  , 


(3) 


■  (4) 


where  is  the  effective  density  of  states  in  the  conduction  band.  Demanding 
that  electrical  neutrality  be  maintained  far  from  the  strain  and  that  V  is  zero 
far  from  the  strain,  i.e.  , 

n  -»  N  and  V  -♦  0  as  x  -*  ®o  , 
then  using  the  relation  between  £  and  p , 

N  =  N  exp  [p/kT]  . 
c 


Therefore,  Equation  (4)  can  be  rewritten  as 
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n  -  N  =  exp  [p/kT]  ^exp  (eV/kT)  -l| 


or. 


n  -  N  =  I  exp  (eV/kT)  -l|  , 


(5) 


where  is  the  carrier  concentration  in  the  absence  of  strain.  Then  combining 
Equations  (3)  and  (5), 


.2..  4iren  . 

^  dx  (eV/kT)  -l|  _  (6) 

For  the  sake  of  simplicity,  S  is  chosen  as 

where  a  is  the  strain  at  the  interface .  Thus,  Equation  (6)  becomes 

_  =  -gsj  +  — ^  I  exp  [eV/kT]  -l|  ,  (7) 

dx  '  ’ 


For  calculation,  the-  potential  is  divided  into  two  parts:  Vj  valid  for  0  <  x  <  a 
and  Vg  valid  for  a  <  x  <  /  .  The  boundary  conditions  are 


^2  = 


V„  and 

2 

dx 

ciV„ 

O 

II 

for 

at  X  =  a, 
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Under  these  conditions,  and  assuming  is  small, 


V 


2 


(8) 


where  w  =  4iren^/e  and  is  an  integration  constant,  hi  the  region  0  <  x  <  a, 
in  order  to  obtain  a  solution  in  closed  form,  it  is  necessary  to  consider  the  strain 
term  dominant.  *  With  this  restriction,  becomes 


= 


+  y  X  +  6 


(9) 


\n4iere  p  =  -{gsa)/a.  and  y  and  6  are  integration  constants.  Evaluating  the 
constants  of  integration,  one  finds 

Vo  =  {(3  a^/2)  [e'^  (1  +  ta)]  , 

-1 

y  =  -/3  a  (3at  +  2)  [2(1  +  ta)] 

6  =  , 
t  =  (we/kT)^^^  . 

3.  Calculation  of  the  Photovoltage 
21  22 

Tauc  ’  has  derived  a  rather  general  expression  for  the  photovoltaic  effect. 
This  expression  is  used  as  a  starting  point  for  the  photovoltage  calculation. 
Tauc  finds  that  the  photovoltage  U  generated  by  an  n-tjiie  semiconductor  is 
given  by : 


'•‘A  consequence  of  this  limitation  is  that  this  model  will  not  be  valid  at  very  high 
light  levels. 
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U 


-W  (  f  Aa  1  dx  +  ~  r  dAg 

e  1  Jo  +  A<r  n  3x  /i  +  /x  Ja  +  A<r 
I  n  p 


(10) 


where  n  is  the  carrier  concentration,  cr  the  conductivity,  Acr  the  change  in  con¬ 
ductivity  with  illumination,  and  Mjj  and  Mp  are  the  electron  and  hole  mobilities, 
respectively.  The  connection  between  the  internal  electrostatic  potential  V  and 
the  quantities  in  Equation  (10)  is  made  by  means  of  the  following  relations 


M  =  ?  -  eV 
d/i  =  0 


n 


N  e 
c 


S/kT 


(11) 


From  these  equations ,  the  relation  between  n  and  V  is  found : 


1  ^  lY 

n  dx  kT  dx 


(12) 


For  convenience,  make  the  following  abbreviatioiiS , 

e  -  e/kT  P  -  dV/dx  , 


and  define  T  (x)  by 


r(x) 


/ 


OPdx 


Then  Equation  (12)  becomes 


(13) 


1  dn  _  dl 
n  dx  ~  dx  ’ 


(13a) 


and 


In  n  =  r(x)  +  constant  , 
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Thus,  for  the  region  a  <  x  <  /,  with  V2  calculated  in  Section  B,  n  becomes 


n  =  n^  exp 


(evykT), 


-mx 


=  n  e 
o 


r2(x) 


(14) 


where  m  =  \  l~  . 


Similarly,  for  the  region  0  <  x  <  a,  n  becomes 


n 


n^  exp 


^e/kT^  (p  x^/2  +  y 


X  + 


Ti  (X) 


n  e 
o 


(15) 


Using  Equations  (13a),  (14)  and  (15),  the  first  term  in  Equation  (10)  can  be 
rewritten 


-kT 


e 


e^  +  Ao/cr 


dr 


where  =  n^  efi^  .  This  integral  is  called  the  chemical  voltage  . 


Substituting  for  and  r2  and  integrating 


Uc  =  ^  'In 
c  e 


1  Au 

hr.  -O _ 

Act 


•+  1 


(16) 


The  second  term  in  Equation  (10)  is  called  the  diffusion  voltage  and  is  some¬ 
times  associated  with  the  Dember  Effect.  This  term  gives  the  contribution  to  the 
photovoltage  due  to  any  difference  in  the  mobilities  of  the  carriers  diffusing  out 
of  the  illuminated  region.  Assuming  that  the  density  of  excess  carriers  falls 
off  linearly  with  distance  outside  the  illuminated  area  one  finds 
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Jsl  ~  \ 


1  + 


M  +  M 


In 


Aa 
_ o 

£111 


Act  I 
1  +  -7-2  \ 

CT(o)  ' 


where  the  illuminated  region  extends  from  0  to  1  and  Act^  is  the  change  in  con¬ 
ductivity  with  illumination  between  0  and  1.  From  Equations  (14)  and  (15) 


CT(0)  and  CT(1)  can  be  found  in  terms  of  Ij  and  F^.  Thus, 


^CT 

_ o 

CT 

O 


-rj(0) 

e 


U.,  =  — 


b-'T 

kT  n  £  , 


1  + 


n  +  n 
n  p 


Act 


1  + 


Combining  Equations  (16)  and  (17),  the  total  photovoltage  U  becomes 


Act 
_ _ o 

CT 

O 

e  Vl  +  l/b/“'  \  Act  -I 


where  b  is  the  ratio  of  electron-to-hole  mobility  and 


ri(0)  = 


As  noted  earlier,  P  =  -(gsQ!)/a  where -oi/a  is  the  slope  of  the  strain.  Thus, 
for  a  linear  slope  (with  point  "a"  fixed) 


n  2  2  dS  . .. 

/3  a  =  gsa  .  (19) 


Equation  (18)  with  (19)  will  be  used  to  express  the  photovoltage  developed  in  a 
single  segment  of  the  crystal.* 

♦Actually  similar  relations  can  be  found  for  any  strain  function  for  which  ^  ^  0 
at  X  =  0  . 
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C.  THE  ADDITION  OF  ELEMENTAL  PHOTOVOLTAGES 

Consider  the  configuration  below,  which  represents  a  ZnS  crystal  composed  of 
three  segments . 


The  two  end  regions  are  identical  and  the  region  b  c  has  slightly  different  pro¬ 
perties  and  the  crystal  is  uniformily  illuminated  from  a  to  d.  The  conduction 
band  edge  in  the  dark  of  such  a  configuration  might  be  represented  by  the  following 
picture . 


In  the  limit  of  very  high  light  levels  it  appears  that  the  voltage  generated  between 
a  and  b  will  be  cancelled  out  by  the  voltage  generated  between  b  and  c,  since  in 
this  limit  the  photovoltage  depends  only  on  the  difference  in  Fermi  level  (chemical 
potential  £  )  between  parts  a  and  c.  In  this  case,  there  can  be  no  addition  of 
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photovoltages  (unless,  of  course,  b  c  is  a  metal).  However,  at  lower  light 

levels,  it  will  be  shown  that  there  can  be  addition  of  photovoltages.  Qualitatively, 

one  notes  from  Equation  (18)  that  the  photovoltage  saturates  at  very  hi^  li^t 

levels.  Thus,  in  order  to  get  addition  it  is  necessary  for  the  photovoltage  in 

22 

b  c  to  be  unsaturated.  Quantitatively,  one  can  express  the  photovoltage  U  as 


U 


2 

1  +  b 


dx 


dx  , 


where  is  the  dark  resistivity  and  is  a  function  of  x.  Since  the  first  and 
third  layers  are  identical,  the  following  relations  hold: 

A(T=Acrj^  a<x<b 

Ao  =  Aff^  b  <  X  <  c 

2 

Ao  =  AOj^  c  <  X  <  d 

For  convenience,  assume  that  b  is  constant.  Then  U  becomes 

^  ^  -  ^o<^>  -  Po<^>]  [  • 

But  from  the  symmetry  of  the  configuration, 


Po(a)  =  "  Po^'^^  ’ 


aiKl  U  becomes 


U 


-kT  2 
e  1  +  b 


2Aorj[p^(b)  -  p^(a)] 
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Thus,  if  AcTj  =  Affg  there  is  no  addition,  but  if  Aff2  Ao^  U  is  nearly 
twice  the  voltage  generated  by  one  segment. 

The  condition  that  Acr  »  Ao^  can  be  met  if  the  absorption  coefficient  of  region 
a  b  is  greater  than  region  b  c.  It  is  maintained  that  at  the  positive  and  negative 
peaks  this  condition  is  fulfilled  for  one  or  the  other  segment,  and  the  individual 
photovoltages  of  alternating  segments  add  to  give  the  observed  high  photovoltages . 
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IV.  EXPERIMENTAL  RESULTS 

A.  DATA  ON  SPECTRAL  RESPONSE 

The  data  to  be  presented  in  the  following  sections  will  justify  the  model  set  up 
in  the  preceding  section.  The  data  will  be  grouped  into  three  categories: 

(1)  Data  which  concern  the  spectral  response  of  the  APE 

(2)  Data  which  concern  the  origin  of  the  voltage  in  individual  layers 

(3)  Other  data  including  results  of  heat  treatment,  y- irradiation,  doping 
and  the  current-voltage  characteristics  of  the  crystals 

Regardless  of  the  mechanism  which  produces  the  voltage  in  a  single  layer  of 
ZnS  (strains,  for  example)  this  model  differs  from  both  Neumark's  and  Tauc's 
models  in  that  it  assumes  that  a  separate  voltage  generator  exists  in  each  seg¬ 
ment.  Tills  enables  the  model  to  account  for  the  anomalously  high  voltages  at 
both  the  positive  and  negative  peaks.  In  the  Neumark  and  Tauc  models,  the  volt¬ 
age  is  generated  between  the  segments,  and  another  mechanism  is  needed  to 
accomplish  the  change  in  polarity  (Tauc's  doping  gradients,  Neumark's  polari¬ 
zation  fields).  Thus,  an  important  part  of  this  model  is  that  the  two  phases  can 
act  independently. 

1.  Spectral  Response 

The  spectral  response  curve  (shown  on  the  following  page)  is  divided  into  six 
regions  for  the  following  discussion,  and  for  convenience  the  longer  wavelength 
polarity  reversal  is  denoted  by  and  the  shorter  wavelength  polarity  rever¬ 

sal  by  . 
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Origin 

Absorption  coefficient  of  hexagonal  region  is 
greater  than  cubic  region . 

Cubic  material  absorption  increases  sharply. 

"Cubic  generators"  have  maximum  output. 

Cubic  material  absorption  decreases  (in  sense 
of  producing  carriers);  hexagonal  absorption  is 
increasing. 

"Hexagonal  generators"  have  maximum  output. 
Hexagonal  material  decreases  in  absorption . 


As  mentioned  in  the  text,  the  transmission  curves  of  hexagonal  and  cubic  mater¬ 
ial  support  the  proposals  for  the  origins  of  Regions  II  and  IV,  that  is,  cubic  ma¬ 
terial  does  begin  to  absorb  in  Region  II  and  hexagonal  material  does  begin  to 
absorb  in  Region  IV. 
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Measurement  of  the  photoconductivity  of  ci^stals  showing  anomalous  effects 
gives  evidence  to  support  the  proposed  origins  of  Regions  IV  and  VI.  The  pro¬ 
posed  decrease  in  carrier  generation  in  the  cubic  phase  in  Region  IV  is  evidenced 
by  the  dip  in  the  photocurrent  (Fig.  10a,  b)  between  335  and  330  m/x .  This  is 
just  to  the  short  wavelength  side  of  the  cubic  band  gap,  and  corresponds  to 
Region  IV.  The  decrease  in  photocurrent  below  310  m/x  supports  the  proposed 
origin  of  Region  VI. 


If  the  origins  of  Regions  I  and  V  are  correct,  then  the  photovoltages  and  currents 
must  be  dependent  on  the  polarization  of  the  incident  light.  It  is  known^^  that 
light  polarized  X  to  the  "c"  axis  is  more  heavily  absorbed  than  light  polarized 
II  to  "c"  (in  hexagonal  material).  Thus,  for  the  X  polarization,  the  voltages 
and  currents  should  be  greater  (in  absolute  value)  than  for  the  ||  polarization. 

This  is  precisely  what  was  reported  by  Cheroff  et  al .  (see  Fig,  11a).  Moreover, 

the  effect  of  polarization  on  the  photovoltage  and  current  in  Region  in  can  also 
be  predicted  from  this  model.  Since  the  hexagonal  "generator''  opposes  the  cu¬ 
bic  "generator, "  then  any  increase  in  the  hexagonal  output  in  this  region  will 
show  up  as  a  decrease  in  the  observed  voltage  in  Region  HI.  Then,  since  cubic 
material  is  isotropic,  the  output  in  Region  in  should  be  greater  for  light  polar¬ 
ized  II  to  "c"  than  for  light  polarized  X  to  "c.  "  This  effect  was  also  observed 
by  Cheroff  et  al.  (Fig.  11a)  and  was  confirmed  by  Beun  and  Goldsmith  (Fig.  11b). 
Figures  11a  and  11b  also  show  a  large  shift  in  ^Qy  'vith  polarization,  while  there 
is  no  shift  in  According  to  this  model,  at  the  cubic  material  is  much 


OR 


more  absorbing  than  the  hexagonal  so  that  no  shift  is  predicted.  On  the  other 
hand,  at  X^y  the  hexagonal  material  is  strongly  absorbing  so  that  one  does  ex¬ 
pect  a  shift  with  polarization. 
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Fig.  10a  Photoconductivity  in  a  Mixed  ZnS  Crystal 


WAVELENGTH  (MILLIMICRONS) 


Relative  photocurrent  of  crystal  D-16.  The  dark  current  at  22.  5  vdc  is  <  0.  001 
X  10“10  amp.  The  photocurrent  at  330  mp  is  ~  0.75  x  10”^®  amp,  giving  a 
relative  change  in  conductivity  of  ~  750.  Arrows  "a"  and  "b"  mark  positions  of 
hexagonal  and  cubic  absorption  edges,  respectively. 
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Fig.  10b  Photoconductivity  in  a  Mixed  ZnS  Crystal 


Relative  photocurrent  of  crystal  D-21.  Arrows  "a"  and  "b"  mark  positions  of 
hexagonal  and  cubic  absorption  edges,  respectively. 
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Fig.  11a  Effect  of  Polarized  Light  on  Photovoltaic  Current 


CHANGE 

OF 

SCALE 


Dependence  of  spectral  response  of  the  photovoltaic  current  on  the  polarization 
of  incident  light,  x  implies  the  electric  vector  E  of  the  incident  light  is  per¬ 
pendicular  to  "c"  axis,  II  implies  E  is  parallel  to  "c"  (after  G.  Cheroff, 

R.  C.  EnckandS.  P.  Keller,  Phys.  Rev.  m.  1091  (1959)], 
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Fig.  lib  Effect  of  Polarized  Light  on  Photovoltaic  Current 


Xifn  fi) 


Dependence  of  spectral  response  of  the  photovoltaic  current  on  the  polarization 
of  incident  light,  x  implies  the  electric  vector  E  of  the  incident  light  is  perpen¬ 
dicular  to  "c"  axis,  II  implies  E  is  parallel  to  "c"  [after  Beun  and  Goldsmith, 
Helv.  Phys.  Acta^,  508  (I960)). 
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2.  Polarity  Reversals 

The  positions  of  the  polarity  reversals  and  are  dependent  on  the  in¬ 
tensity  of  the  incident  light  (Fig.  12a  and  b  are  typical) .  hi  fact,  in  some  crys¬ 
tals  the  polarity  reversals  dis^pear  altogether  at  low  light  levels .  The  data 
presented  in  Fig.  12a,  b  were  taken  by  adjusting  the  monochromator  and  light 
source  for  maximum  short  circuit  current  at  a  convenient  wavelength  and  then 
decreasing  the  light  level  by  placing  copper  screens  in  the  light  beam.  In  order 
to  analyze  Fig.  12a,  b  in  terms  of  the  proposed  model,  it  is  necessary  to  know 
quantitatively  the  variation  of  the  absorption  constants  with  wavelength  in  this 
region.  However,  this  same  phenomenon  can  be  seen  in  another  experiment. 
For  Fig.  13a,  b,  the  voltages  at  ^Qy  and  ^•qj^  were  measured  as  the  intensity 
was  decreased.  Of  course,  the  increase  in  voltage  noted  does  not  continue  in¬ 
definitely,  since  the  photovoltaic  effect  must  disappear  for  zero  light  intensity. 

This  variation  in  voltage  at  and  X^^  can  be  analyzed  in  terms  of  the  two- 
segment  unit  model.  Let  the  voltage  of  a  two-segment  section  of  the  crystal  be 
represented  by  the  following  equation 


U 

n 


(  1  +  i  (  1  +  ( 

""ll  +  rjiluij  -  l"  |H-r,2lu2^ 


(20) 


where  is  expressed  in  units  of  kT/e,  I  is  the  maximum  light  intensity  in  the 
experimental  set-up,  r]  and  ri  are  constants  of  the  cubic  and  hexagonal  seg- 
ments,  respectively.  Vl  is  defined  by  the  relation 


ijl 


a 

o 


Uj^  and  Ug  are  related  to  the  fields  in  the  crystal  due  to  the  stra.'ns  and  are  in¬ 
dependent  of  light  intensity.  At  ^Qy  »  for  example, 
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Fig.  12a  Dependence  of  Red  Polarity  Reversal  on  Light  Intensity 
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RELATIVE  LIGHT  LEVEL  I/Io  (ARBITRARY  UNITS) 


Dependence  of  the  wavelength  of  red  polarity  reversal  on  the  intensity  of  incident 
light.  The  maximum  light  intensity  available  in  the  experiment  is  arbitrarily  de¬ 
noted  1.  0. 
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Fig.  12b  Dependence  of  Violet  Polarity  Reversal  on  Light  Intensity 


RELATIVE  LIGHT  LEVEL  I/I^  (ARBITRARY  UNITS) 


Dependence  of  the  wavelength  of  violet  polarity  reversal  on  the  intensity  of  inci¬ 
dent  light.  The  maximum  light  intensity  available  in  the  experiment  is  arbitrar¬ 
ily  denoted  1. 0. 


LOCKHEED  MISSILES  &  SPACE  COMPANY 


OPEN  CIRCUIT  VOLTAGE  AT  XqVi  tVOLTS) 


44 


Fig.  13a  Experiaiental  and  Theoretical  Dependence  of  at  on 

Light  Intensity 


Experimental  and  theoretical  dependence  of  Vqc  ^OV  on  intensity  of  inci¬ 
dent  light.  Open  circles  are  experimental  points;  filled  circles  are  calculated. 

Calculated  from: 

^OC  "  ^”G.4  -  0.55)  ’ 
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Fig.  13b  Experimental  and  Theoretical  Dependence  of  Vqc 
at  XqV  on  Light  Intensity 


Experimental  and  theoretical  dependence  of  Vqc  at  Xqv  on 
intensity  of  incident  light.  Open  circles  are  experimental 
points,  and  filled  circles  are  calculated. 

Calculated  from; 

'^oc  ^  ^"(i  -  o.se)  ’ 
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Therefore , 


Un  =  0 


In 


1  +T)^I 

1  +  T)  jUjI 


l+r,2l 


1  +  JJj^I  1  +  ’^2^ 


(21) 


Solving  this  relation  for  n  „  , 


V 


2 


r»il 

- 

1 

t?il| 

("2 

-“0 

1  + 

“2  -  ^ 

at  Xq  .  The  following  part  of  the  derivation  is  valid  for  X^y.  In  this  region 
Aa/Oo  for  the  cubic  phase  has  decreased  so  that  0^1  is  small.  Moreover  both 
U2  and  Uj  are  of  the  order  of  0.  01,  *  so  that  -  u^)  is  negligible  compared 

to  U2  -  1  .  Therefore, 


ri 


2  * 


M  ^ 

Uo  -  1 


(22) 


Now,  decrease  the  light  level  to  1-6  and  recalculate  at  X^^  using  Equation 
(21).  Calling  this  voltage 


(1  +  ’Ji(I-6)+17jT)2U2(I-6)^ 

|l  +  ■n^{l  -  6)  +Tj^u^(I  -6)  +  tj^172Uj^(I  -  6)^ 


(23) 


Then,  neglecting  terms  in  6 


2 


and  using  relation  (21),  becomes 


♦The  voltage  generated  per  layer  should  be  about  0. 1  volts'^  or  4  (kT/e)  at  room 
temperature.  Thus  (maximum)  =  4  or  lnu-1  =  4.  Thus  u  =  e"4  s  o.Ol. 
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1  +  ^2^  +  TiiUiI  +  -  6(r,i  +  -  Zy^n^u^l)  ] 

I  1  +  ijgl  ■*■  "  *^(^2  '*’  ^l'*!  "  ^ 


The  equation  has  the  form 


U’ 

n 


In 


(a  -  b6  ( 
j  a  -  c6  j 


where 

b  =  ^,^  +  ^2^2  -2^i’^2V  I 

c  =  7,2  +  Tj^u^  -  arj^Tj^u^I  j  . 
Then,  using  Equation  (22)  and  Equations  (25), 


(24) 


(25) 


(26) 


but,  following  the  reasoning  leading  to  Equation  (22),  Equation  (26)  can  be  re¬ 
written  as 


-^l)  _  ^  -  ^1  -  ^2 


This  equation  is  symmetric  in  u^  and  u^  and  since  both  u^^  and  Ug  are  <  1, 
it  is  easily  seen  that  ^  >  1. 

Thus, 


0  < 


a  -  bd 
a  -  cd 


<  1 
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for  small  6  ,  and 


i . e the  voltage  goes  negative  as  5  increases. 

The  best  fitting  expressions  of  the  form  of  Equation  (24)  are  shown  on  Fig.  13a, b. 
These  are  just  empirical  fits  and  little  significance  can  be  tied  to  the  various 
values  of  a,  b  and  c  used  in  these  equations.  It  must  be  remembered  that  in 
this  derivation,  only  the  voltage  due  to  two  segments  of  the  crystal  was  calculated. 
Actually,  the  crystal  is  made  of  many  such  layers,  all  slightly  different.  Thus, 
the  experimental  data  is  just  the  resultant  of  the  actions  of  many  slightly  differ¬ 
ent  sets  of  segments . 

3 .  Dependence  of  X  on  the  Elemental  Photovoltage 

According  to  the  proposed  model,  one  would  expect  that  the  relative  magnitudes 
of  the  positive  and  negative  peaks  would  have  an  effect  on  the  position  of  ^Qy 
Qualitatively,  if  the  positive  peak  is  much  larger  than  the  negative  peak,  then 
the  positive  peak  voltage  should  persist  to  shorter  wavelengths.  Such  a  crystal 
would  have  X  ^y  much  farther  in  the  violet  than  a  crystal  in  which  the  positive 
and  negative  peaks  are  of  comparable  magnitude.  The  wavelength  ^Qy  versus 
the  ratio  of  the  magnitudes  of  the  positive  and  negative  peak  voltages  is  shown 
in  Fig.  14  for  a  number  of  crystals.  The  scatter  of  the  data  is  large,  but  one 
observes  a  definite  trend  of  the  type  outlined  above . 

This  trend  has  a  quantitative  basis.  To  show  this  relation,  assume  that  at  the 
positive  and  negative  peaks  the  observed  voltage  is  saturated.  Then,  the  observed 


LOCKHEED  MISSILES  &  SPACE  COMPANY 


6-90-63-20 


49 


positive  and  negative  peak  voltages  and  V_  ,  respectively,  are 


V  =n!^ln-i  . 
+  e  Uj 


„  kT  ,  1 

V  =  -  n  —  In  —  , 
e  U2 


where  n  is  just  half  the  number  of  segments  in  the  crystal,  and  u^^  and  u^ 
are  functions  depending  only  on  the  strain  and  piezoelectric  constants .  Also, 
let  r  be  defined  as  the  ratio  of  the  magnitudes  of  and  V 

In  1/u^ 

’  |v:|  "  CHTJ;  ■ 


then 


r 


and 


=  Ug  .  (28) 

But,  since  we  are  interested  only  in  Qy  .  use  condition  (22)  from  the  preceding 
section. 


Then, 


(29) 
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2 

since  =  j-z —  ,  ^7  o  is  proportional  to  the  absorption  coefficient  of  hexa- 
^  *^20 

gonal  ZnS  and  varies  inversely  with  wavelength.  In  fact  Beun  and  Goldsmith^^ 

2 

find  that  the  absorption  coefficient  a  varies  with  (hv)  .  Thusi 


V 


2 


1 

\ 


2 


and  Equation  (29)  becomes 


X 


2^ 


or, 


r  being  always  greater  than  zero.  Since  r  =  0  implies  =  0  ,  this  means 
that  there  is  no  sign  reversal  and  x  is  not  defined.  Figure  15  is  a  plot  of 
Equation  (30)  for  various  values  of  u„  and  should  be  similar  in  form  to  the  data 
presented  in  Fig.  14.  hi  Fig.  14,  u^  is  not  constant.  Figure  15  shows  that  the 
largest  spread  between  curves  with  different  u^  occurs  for  small  values  of  r. 

This  is  ^parently  the  case  in  Fig.  14.  Also  the  trend  predicted  by  Equation  (30) 
shows  up  in  Fig.  15.  Admittedly,  the  agreement  is  not  excellent,  but  it  must 
be  remembered  that  the  calculation  was  done  for  an  ideal  model  consisting  of 
n  identical  segments  of  hexagonal  and  n  identical  segments  of  cubic  material. 

Visual  observation  of  ZnS  crystals  reveals  that  in  a  real  crystal  the  sizes  of  the  layers 
are  rather  random.  Thus,  the  scatter  in  points  on  Fig.  14  is  hardly  unexpected. 
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I 

I 

I. 

! 

I  Fig.  14  Relation  Between  and  r 


i  The  wavelength  of  the  violet  polarity  reversal  vs.  the  absolute  magnitude  of  the 

*  ratio  of  the  maximum  open-circuit  voltage  at  the  positive  peak  to  the  maximum 

value  of  the  open-circuit  voltage  at  the  negative  peak.  This  ratio  is  denoted  r. 

I 
1 
1 

I 
I 
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Furthermore,  Equation  (30)  was  calculated  assuming  that  the  positive  and  negative 
peak  voltages  are  saturated.  Experimentally,  it  was  not  always  possible  to 
saturate  the  photovoltages. 

4.  Temperature  Dependence  of  the  Spectral  Response 
28 

Cheroff  has  measured  the  temperature  dependence  of  the  spectral  response. 

Two  pieces  of  data  are  forthcoming.  First,  that  the  voltage  Increases  at  low 
temperatures  according  to  the  relation 


V 


OC 


Second,  the  difference  in  energy  AE  between  and  X^j^  increases  linearly  with 

temperature.  The  slope  at  low  temperature  is  AEAT  =  3.7,  and  above  room 

temperature  the  slope  is  2AE/kT.  According  to  the  proposed  model,  Xqj^  is 
associated  with  the  cubic  band  gap  and  X  is  associated  with  the  hexagonal 
band  gap.  Thus,  the  change  with  temperature  of 

^  _  he  he 
^OR  ^OV 

is  just  given  by  the  difference  in  d  Eg/dT  between  the  cubic  and  hexagonal 
phases  and  should  correspond  to  the  slope  of  the  AE  vs.  T  curve  of  Cheroff. 

Above  room  temperature  this  difference  is  just  AE  =  2kT  (Ref.  29). 

B.  EVIDENCE  OF  ORIGIN  OF  PHOTOVGLTAGE 

Thus  far  the  evidence  presented  has  only  established  the  fact  that  the  spectral 
response  of  the  photovoltage  can  be  accounted  for  on  the  basis  of  the  proposed 
model.  However,  nothing  has  been  said  of  the  origin  of  the  voltage  in  each 
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segment.  The  following  sections  will  present  evidence  that  the  elemental 
photovoltages  are  due  to  piezoelectric  fields  induced  by  permanent  non-uniform 
internal  strains. 

1.  Measurements  on  Mixed  Crystals 

There  is  some  very  qualitative  evidence  that  the  voltage  generators  are  associ¬ 
ated  with  some  property  of  the  crystal  structure.  It  is  found  that  in  the  majority 
of  crystals  studied,  the  end  of  the  crystal  grown  last  is  positive  when  illuminated 

4 

with  light  of  wavelength  ~340mM  .  This  result  is  supported  by  Merz.  More¬ 
over,  it  was  also  found  that  when  a  compressional  strain  is  applied  in  the  dark 
along  the  long  axis  of  the  crystal,  a  piezoelectric  voltage  is  generated.  The 

end  of  the  crystal  which  becomes  positive  under  compression  is  identical  with 

17 

the  end  of  the  crystal  which  is  positive  at  340  m/i.  Lempicki  found  a  one-to- 
one  correspondence  between  the  polarity  of  the  positive  peak  and  the  pyroelectric 
voltage. 

In  order  to  test  the  postulate  that  strain-induced  piezoelectric  fields  are  the 
voltage-producing  mechanisms  in  the  crystal,  the  spectral  response  of  the  APE 
was  measured  on  crystals  with  a  longitudinally  applied  stress.  The  stress  was 
always  compressional,  since  the  crystals  are  too  weak  to  stress  in  tension. 

To  analyze  the  effect  of  application  of  an  applied  stress  on  the  crystal,  consider 
the  figure  on  the  following  page.  The  solid  line  represents  the  variation  of  the 
transverse  lattice  constant  on  passing  from  one  phase  to  the  other.  Suppose  a 
uniform  stress  IT  is  applied  to  the  crystal.  Then  in  the  bulk  of  each  segment,  the 

lattice  dilates.  The  dilation  in  Region  I  is  and  in  Region  n  is  5^ 

assume  a  linear  gradation  In  lattice  constant  from  the  hexagonal  to  cubic  phase). 
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Then,  if  this  slope  before  application  is 

Aa/A/ 

then,  after  an  applied  stress  it  is 


Aa' 


A/ 


But,  since  the  stress  is  uniform. 


where  c-  and  c„  are  the  appropriate  elastic  constants.  Thus. 
1  ^ 


,  Aa  +  6^  -  -1 
^  =  \  2  / 
A  A/ 
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The  stress  11  is 


and 


therefore. 


n 


Vi 


6j  =  n/cj^ 


Aa 

A/ 


Aa.n^,-^ 

A/ 


Unfortunately,  values  for  both  and  are  not  available.  Thus,  it  cannot 
be  predicted  whether  the  slope  increases  or  decreases. 

With  that  difficulty  in  mind,  some  analysis  of  the  effect  of  strain  can  be  made. 
As  in  the  section  on  intensity,  consider  a  two-segment  section  of  the  crystal 
and  represent  the  resultant  voltage  as 


U  = 
n 


1  +  a, 


In 


-  <y. 


1  +  Sj^e 


1  +  a2C 
1  +  a„ 


where  U  is  in  units  of  kT/e  ,  a-  and  a  are  Ao/cr  for  each  phase,  and 
n  _  1  o 

-CL  - 

e  and  e  are  the  functions  depending  on  strain  gradient,  etc.  Then 
apply  a  small  stress  p  .  becomes 


U  =  In 
n 


1  +  a 


-C2+ 


where  and  are  constants. 
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At  U  -  0  so  that 

OV’  n 


-ev 


1  +  ai  /  1  +  a2e 


^  *  “2  \  ,  -‘'l 

1  +  a^e 


=  1 


Then,  if  the  applied  stress  is  small, 


e  =  1  +  , 


and  can  be  rewritten 


Un  =  in 


+  a2e 


^(l  +  a2)(l  +  a^e  ^)  +  a^e 


(31) 


This  equation  has  the  form 


i„  (4-1^ 

\  A  +  cp 


(32) 


The  argument  of  Expression  (32)  has  the  form 
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Thus,  for  either  case  will  saturate.  This  is  just  what  was  found.  The  best 
data  taken  at  is  shown  in  Figs.  16  and  17.  As  with  the  light  intensity 

data,  the  results  are  presented  in  two  ways:  (1)  as  a  shift  in  X  Fig.  16, 
and  (2)  as  a  change  in  voltage  at  ^Qyi  Fig.  17. 

2.  Experiments  on  Single  Crystal  ZnS 

Since  the  results  of  the  last  section  show  an  uncertainty  in  the  predicted  polarity 
of  the  voltage  generated  at  by  an  applied  uniaxial  stress,  an  experiment 
was  devised  using  single  crystal  ZnS  to  test  directly  the  origin  of  the  photovolt¬ 
age  in  one  layer  of  a  mixed  crystal.  The  experiment  was  preformed  on  a  single 
crystal  of  hexagonal  ZnS.  A  similar  experiment  was  attempted  on  cubic  ma¬ 
terial,  but  it  was  found  that  the  resistivity  of  the  material  available  was  so  low 
that  the  crystal  showed  no  piezoelectric  properties. 

The  experiment  consisted  of  non-uniformly  stressing  the  crystal  and  then  ob¬ 
serving  the  photovoltage  generated  when  the  strained  crystal  was  illuminated 
with  ultraviolet  li^t.  The  contribution  to  the  photovoltage  due  to  contacts  was 
taken  into  account  by  noting  the  photovoltage  in  the  unstrained  crystal  and  sub¬ 
tracting  off  this  effect. 

In  order  to  interpret  the  results  of  this  experiment,  one  needs  a  model  of  the 
strain  pattern  in  the  mixed  crystal.  This  subject  was  touched  on  in  Section 
in-B,  where  it  was  pointed  out  that  the  hexagonal -to-cubic  transition  is  much 
less  abrupt  than  the  cubic-to-hexagonal  transition,  and  that  the  strain  gradients 
will  be  smaller  at  the  hexagonal -to-cubic  transition.  It  will  be  assumed,  in 
keeping  with  the  qualitative  nature  of  this  work,  that  these  photovoltages  are 
negligible  in  comparison  to  the  photovoltages  generated  at  the  cubic-to-hexagonal 
transition.  Thus,  in  a  mixed  crystal,  the  photovoltages  associated  with  the 
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Fig.  16  Shift  of  With  Applied  Stress 


APPLIED  STRESS  (GRAMS) 


The  shift  in  position  of  the  violet  polarity  reversal  with  compressional  stress 
applied  along  the  "c"  direction. 
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Fig.  17  Open-Circuit  Voltage  at  as  a  Function  of  Applied  Stress 


20  40  60  80  100 

APn.lED  STRESS  (GRAMS) 


The  change  in  voltage  at  the  violet  polarity  reversal  with  compressional  stress 
applied  along  the  "c"  direction. 
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hexagonal  layers  will  be  due  to  a  strain  gradient  which,  in  each  hexagonal 
layer,  decreases  in  the  direction  of  crystal  growth. 

4 

It  was  pointed  out  in  Section  I-C  that  Merz  found  a  relation  between  the  direc¬ 
tion  of  crystal  growth  and  the  polarity  of  the  photovoltage  at  340  ran.  This  re¬ 
sult  was  confirmed  for  the  vast  majority  of  crystals  examined  in  this  study  and 
it  was  found  that  the  end  of  the  crystal  grown  last  was  always  positive  when 

illuminated  with  light  of  wavelength  corresponding  to  the  positive  peak. 

17 

Lempicki's  result  that  the  end  of  the  crystal  which  became  positive  on  heating 
(secondary  pyroelectric  effect)  was  also  positive  when  illuminated  with  light  of 
wavelength  340  nut  (the  positive  peak)  is  also  necessary  to  the  interpretation 
of  this  e;q)eriment.  Since  only  the  hexagonal  phase  is  pyroelectric,  this  ef¬ 
fect  allows  one  to  orient  the  hexagonal  single  crystal,  relative  to  its  orientation 
in  the  mixed  crystal. 

To  summarize:  it  is  now  evident  that  the  hexagonal  layers  in  the  mixed  crystals 
are  so  oriented  that  the  strain  gradients  in  these  layers  decrease  toward  that 
end  of  the  crystal  which  becomes  positive  when  illuminated  with  light  corres¬ 
ponding  to  the  positive  peak.  Furthermore,  the  "c"  axis  is  so  oriented  that  on 
heating  (secondary  pyroelectric  effect)  the  end  of  the  crystal  which  becomes 
positive  is  identical  with  the  end  which  becomes  positive  when  the  crystal  is 
illuminated  with  li^t  corresponding  to  the  positive  peak.  Thus,  in  the  single 
crystal,  the  strain  should  decrease  in  the  direction  of  the  secondary  pyroelec¬ 
tric  effect  on  heating. 

In  carrying  out  the  experiment,  the  direction  of  the  secondary  pyroelectric  ef- 

31 

feet  was  found  by  using  the  dynamic  method  described  by  Chynoweth.  When 
this  direction  was  identified,  the  crystal  was  ground  to  the  shape  shown  on  the 
following  page. 
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DIRECTION  OF 
SECONDARY 
PYROELECTRIC 
EFFECT  ON 
HEATING 


DIRECTION  OF 

DECREASING 

STRAIN 


A  stress  was  then  applied  perpendicular  to  the  "c"  axis  and  perpendicular  to  the 
trapezoidal  face  of  the  crystal.  In  this  manner,  a  strain  gradient  of  the  trans¬ 
verse  lattice  constants  can  exist  aloi^  the  "c"  axis  when  the  center-of-mass  of 

the  load  is  not  directly  above  the  centroid  of  the  trapezoidal  face  of  the  sample. 

32 

For  this  case,  the  stress  gradient  da/dx  is  given  by  the  relation 

da  F  g 
dx  "  I 

XX 

where  F  is  the  applied  force,  i  the  distance  on  the  face  of  trapezoid  between 
the  center-of-mass  of  the  load  (b)  and  the  centroid  of  the  sample  (a).  is 
the  moment  of  inertia  about  the  line  x-x.  The  photovoltage  was  measured  be¬ 
tween  faces  d,  and  d„  .  This  configuration  was  used  instead  of  one  in  which 
the  stress  was  applied  along  "c,"  so  that  the  observed  photovoltage  would  involve 
only  the  transverse  piezoelectric  constants  and  not  the  longitudinal  piezoelectric 
constants . 

When  the  stress  was  applied  and  the  sample  illuminated,  a  net  pholovoltage  of 
^proximately  4  mv  was  observed.  The  polarity  was  such  that  contact  d^  was 
positive  with  respect  to  62  ■  That  is,  the  polarity  is  opposite  the  polarity  as¬ 
sociated  with  the  secondary  pyroelectric  effect  on  heating,  and  therefore,  the 
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polarity  associated  with  the  positive  peak  photovohage.  This  is  precisely  the 
required  result,  since,  in  the  model,  the  photovoltage  due  to  the  hexagonal 
layers  in  the  mixed  crystal  is  associated  with  the  negative  peak. 

The  value  of  the  photovoltage  one  expects  in  this  experiment  can  be  calculated 
using  the  expression  derived  for  U  in  Section  ni-B.  2.c.  From  this  equation 
the  maximum  value  of  the  photovoltage  is  found  to  be 

=  gsa^  d  n/dx  . 

The  piezoelectric  constant  g  was  measured  on  the  same  sample  used  in  this 

_9 

e:qx!riment  and  was  found  to  be  2. 7  x  10  (v-cm/dyne).  *  Using  this  value  for 
g  the  maximum  value  of  the  photovoltage  expected  for  the  configuration  used 
in  the  experiment  is 


U  =  3.4  X  10"^  volts  . 
m 

This  value  appears  to  be  slightly  too  low  but  it  must  be  remembered  that  a 

strain  gradient  can  also  produce  a  photovoltage  by  virtue  of  the  band  gap  change 

22 

produced  by  the  strain  gradient.  This  effect  was  investigated  by  Tauc.  On 
the  basis  of  his  calculation  and  using  a  value  of  dEg/dP  from  Ref.  29,  this 

_3 

contribution  to  the  photovoltage  is  -0.6  x  10  volts.  Note  this  voltage  is  of  the 

opposite  polarity  from  the  voltage  observed.  Thus,  including  this  correction 

in  the  calculation  of  the  photovoltage  in  a  strained  crystal,  one  finds  that  the 

-3 

observed  photovoltage  should  be  4 x  10  volts.  This  is  just  the  voltage  which 
was  observed. 

’Tor  the  sake  of  completeness,  the  longitudinal  piezoelectric  constant  g„-  was 
measured  and  found  to  be  5. 1  x  10'®  (v-cm/dyne). 
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The  agreement  between  the  two  voltages  is  probably  better  than  is  warranted, 
due  to  errors  involved  in  the  experiment.  For  instance,  the  quantity  i  in  the 
stress  gradient  expression  was  very  difficult  to  measure  and  could  be  in  error 
by  50%.  Thus,  the  agreement  is  probably  no  better  than  order-of-magnitude 
agreement.  However,  the  experiment  still  shows  that  the  observed  polarity  of 
the  photovoltage  is  of  the  right  polarity  to  account  for  the  voltage  at  the  negative 
peak  and  that  the  calculations  in  Section  III-B.  3  give  better  than  order-of- 
magnitude  agreement  with  the  experiment. 
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V.  OTHER  RESULTS 

This  section  will  be  devoted  to  the  presentation  and  discussion  of  experiments 
which  yield  results  that  were  not  reproducible.  That  is,  either  the  experiment 
was  not  repeatable  on  the  same  crystal  ax>d/or  with  different  crystals  (with 
similar  photovoltaic  properties).  Since  the  magnitude  of  the  effects  observed 
varied  greatly  between  crystals,  it  was  felt  necessary  that  measurements  should 
always  be  repeated  for  a  number  of  crystals.  By  and  large,  the  data  presented  in 
this  section  are  the  results  of  experiments  in  which  no  consistence  could  be  foimd 
between  crystals.  The  e^qieriments  which  fell  into  this  category  were  as  follows; 

1)  The  Study  of  Surface  Conditions 

2)  The  Effect  of  y -radiation  on  the  APE 

3)  Effects  of  Diffusion  of  Ag  into  Crystals  and  Heat  Treatment 

4)  Effects  of  Addition  of  Cu  to  Starting  ZnS 

5)  Current  Voltage  Characteristics  of  Crystals  Showing  APE 

A.  THE  STUDY  OF  SURFACE  CONDITIONS 

In  an  attempt  to  study  the  role  played  by  surface  conditions  on  the  APE,  a 
number  of  crystals  whose  spectral  response  characteristics  were  known  were 
subjected  to  a  mild  sandblasting  perpendicular  to  the  "c"  axis,  using  micron¬ 
sized  aluminum  oxide  powder.  After  the  first  treatment,  most  of  the  crystals 
showed  a  decrease  in  the  short-circuit  current.  The  open-circuit  voltage  either 
remained  constant  or  increased.  The  overall  shape  of  the  spectral  response 
curve  appeared  to  be  unchanged.  The  decrease  in  the  short-circuit  current  is 
attributable  to  a  loss  in  the  number  of  hole -electron  pairs  created  by  the  in¬ 
cident  light.  The  decrease  in  the  number  of  hole -electron  pairs  may  be  attri¬ 
buted  to  either  one  or  both  of  the  following  causes.  First,  the  reflection  losses 
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at  the  surface  will  be  increased  by  the  irregularities  produced  on  the  surface 

by  the  grit.  Second,  since  the  surface-to-volume  ratio  of  a  typical  crystal 

is  about  50  cm~^,  surface  recombination  might  be  expected  to  be  important  in 

determining  the  lifetime  of  excess  carriers.  Generally,  sandblasting  a  surface 

increases  the  surface  recombination  and  thus,  decreases  the  lifetime  of  the 

excess  carriers.  Since  in  the  steady  state  the  number  of  excess  carriers  N 

33 

is  related  to  the  lifetime  r  and  the  generation  rate  f  by 


N  =  fT, 


a  decrease  in  the  excess  carrier  concentration  will  result  from  a  decrease  in 
T  .  This  leads  to  a  decrease  In  the  short-circuit  current 

The  increase  in  the  open-circuit  voltage  could  be  due  to  an  Increase  in  the 
amount  of  disorder  in  the  crystal  (i.e. ,  the  number  of  segments).  It  has  been 
found  that  even  mild  pressure,  ^plied  at  a  point  perpendicular  to  the  "c"  axis. 
Induces  additional  disorder  in  the  vicinity  of  the  stress.  *  In  fact,  a  number  of 
crystals  developed  cracks  transverse  to  the  "c"  axis  as  a  result  of  the  sand¬ 
blasting.  There  was  no  consistency  in  the  results  of  further  sandblasting, 
except  that  most  crystals  developed  transverse  cracks. 

B.  THE  EFFECT  OF  y  -  RADIATION  ON  THE  APE 

60 

A  number  of  crystals  were  irradiated  with  Co  gamma  rays .  Two  irradiations 

5  4 

were  made,  one  at  a  flux  density  of  5x  10  r/hour  and  the  other  at  5x  10  r/hour. 

A  e 

The  integrated  doses  were  10°  r  and  10  r,  respectively.  The  purpose  of  these 
irradiations  was  to  try  to  study  the  effect  of  impurities  on  the  APE. 


*J.  Singer,  Private  Communication,  General  Telephone  Co. 
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60 

Co  gamma  rays  will  produce  vacancy- interstitial  defects  which  can  act  as  donors 

g 

or  acceptors .  After  irradiation  to  10  r,  the  crystals  displayed  no  anomalous 
photovoltaic  effects.  That  is,  high  photovoltages  and  polarity  changes  were 
absent.  In  addition,  it  was  noted  that  the  conductivity  of  these  crystals  increased 

5 

by  factors  of  the  order  of  10  .  This  increase  in  conductivity  could  well  be  the 

cause  of  the  disappearance  of  the  APE.  As  the  model  is  presented  here, 

this  connection  is  not  evident.  This  follows  from  the  approximations  made  in 

the  derivation  of  the  photovoltage  for  a  single  segment  of  the  crystal.  It  was 

necessary  to  makf'  ♦^’ese  approximations  in  order  to  obtain  a  solution  in  closed 

form.  However,  it  is  well  known  that  a  crystal  of  low  resistivity  will  show  no 

34 

piezoelectric  properties.  Thus,  the  disappearance  of  the  APE  in  these 
crystals  may  well  be  due  to  the  shorting  out  of  the  piezoelectric  properties  of 
these  crystals  by  the  increased  conductivity. 

35 

A  calculation  of  the  Debye  length  for  irradiated  and  unirradiated  crystals 

gives  some  substance  to  this  argument.  For  a  typical  crystal (p  ~  10^^  S2  -  cm), 

the  Debye  length  is  approximately  0.2  cm.  However,  in  an  irradiated  crystal, 

the  Debye  length  is  decreased  by  a  factor  of  about  (  to  a  value  of  ap- 

-4  4 

proximately  7x10  cm.  Merz  estimates  that  the  segments  in  an  anomalous 

-3 

crystal  are  of  the  order  of  10  cm  long.  Thus,  the  Debye  length  in  an  irradi¬ 
ated  crystal  is  of  the  order  of,  or  less  than,  the  length  of  a  representative 
segment  of  the  crystal.  This  makes  the  presence  of  a  piezoelectric  field  in 
such  a  segment  somewhat  unlikely. 

5 

Among  the  crystals  which  were  irradiated  to  10  r,  the  effects  were  not  as 
drastic.  In  all  cases,  the  short-circuit  current  increased  in  the  vicinity  of  the 
negative  peak.  In  both  the  positive  peak  region  and  the  long  wavelength  region 
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(X  >  X^„),  the  behavior  was  varied.  There  appeared  to  be  no  correlation 
between  changes  in  current  in  one  region  with  changes  in  another. 

The  open-circuit  voltage  also  changed  with  irradiation.  It  was  found  that  if  the 
voltage  at  wavelengths  corresponding  to  the  negative  peak  increased,  the  volt¬ 
age  at  the  positive  peak  decreased.  The  converse  was  also  true.  This  con¬ 
forms  to  the  model,  in  that  it  is  postulated  that  the  generators  in  adjacent 
layers  are  opposed.  Thus,  for  the  purposes  of  discussion,  suppose  that  the 
generators  in  the  hexagonal  regions  increase  their  output.  Then  in  the  region 
of  the  positive  peak  (i.e. ,  the  region  assoeiated  with  the  eubic  segments),  the 
photovoltage  should  decrease,  since  any  contribution  from  the  hexagonal  layers 
is  of  the  opposite  polarity.  Furthermore,  the  absolute  magnitude  of  the  change 
in  voltage  at  the  positive  peak  is  alw'ays  close  to  one-half  the  change  in  voltage 
at  the  negative  peak. 

The  changes  in  the  positions  of  the  polarity  reversals  were  not  consistent. 

C.  EFFECTS  OF  DIFFUSION  OF  Ag  INTO  CRYSTALS  AND  HEAT  TREATMENT 

In  another  attempt  to  study  the  effects  of  impurities  on  the  APE,  silver  (in 

the  form  of  silver  sulfide)  was  diffused  into  crystals  with  known  spectral  re- 
^  fi 

sponses.  The  diffusion  was  carried  out  at  about  600 °C.  The  results  were 
entirely  erratic,  most  of  the  crystals  losing  their  anomalous  properties.  The 
conditions  of  the  experiment  were  not  well  controlled,  as  it  was  necessary  to 
remove  the  electrical  contacts  from  the  samples  before  the  diffusion.  Since  the 
crystals  are  small,  it  was  not  possible  to  replace  the  contacts  in  exactly  the  same 
place.  This,  in  itself,  can  change  the  spectral  response  of  the  crystals  (since 
different  segments  may  be  shorted  out  by  the  electrode  material).  Furthermore, 
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it  is  known  that  heat  treatment  can  change  the  structure  of  ZnS  (Ref.  37).  Thus,  it 
was  not  possible  to  separate  the  effects  of  heat  treatment  from  the  effects  of  doping 


D.  EFFECTS  OF  ADDITION  OF  Cu  TO  STARTING  ZnS 

A  number  of  crystals  were  grown  from  ZnS  to  which  Cu  had  been  added.® 

Several  batches  of  crystals  were  grown  with  doping  concentrations  ranging  from 
-3  -5 

10  to  5x  10  gm  Cu  per  gm  of  ZnS.  Data  was  taken  on  approximately  ten 
crystals  from  each  batch.  In  general,  the  spectral  response  curves  of  crystals 
from  different  batches  were  very  similar.  No  correlation  was  found  between 
the  impurity  concentration  in  the  starting  ZnS,  any  significant  experimental 
quantity  (such  as  ^Qy.  ^qr^’  position  and  magnitude  of  the  positive 

and  negative  peaks .  It  Ciumot  be  stated  unequivocally  that  the  concentration  of 
impurities  in  the  crj’stal  s  was  even  proportional  to  the  concentration  of  im¬ 
purities  in  the  starting  niaterial.  On  this  basis  it  is  felt  that  no  conclusions 
can  be  drawn  from  this  data. 


E.  CURRENT- VOLTAGE  CHARACTERISTICS  OF  CRYSTALS  SHOWING  APE 

The  current-voltage  characteristics  of  crystals  having  similar  photovoltaic 
properties  fell  into  two  classes.  None  of  the  current -voltage  characteristics 
depended  on  the  material  used  to  make  electrical  contact  to  the  crystals.  The 
applied  voltages  ranged  from  millivolts  to  280  volts.  For  very  low  voltages 
(  <  10  volts)  most  crystals  displayed  ohmic  characteristics .  At  higher  voltages 
the  behavior  is  varied.  The  two  categories  of  behavior  are  as  follows: 

1)  Ohmic  to  dependence 

2)  Current  saturation  at  high  voltage 
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1 .  Square  Law  Behavior 

Only  two  samples  displayed  current-voltage  (I  -  V)  characteristics  expressible 
by  an  equation  of  the  form 


I  =  (constant)  v”  . 

The  spectral  response  of  the  APE  in  these  crystals  was  of  the  type  shown  in 
Fig.  1.  One  sample  (L-10)  remained  ohmic  to  40  volts  (^proximately  200 
v/cm).  Above  40  volts  the  current  varied  as  V^.  The  other  sample  (J-5)  re¬ 
mained  ohmic  to  about  the  same  voltage  and  then  changed  to  a  behavior  expres- 
sible  as  I  =  (const.)  V  The  current -voltage  characteristics  of  sample 

(L-10)  is  that  of  a  trap-free  solid  whose  high  field  conductivity  is  limited  by 
space  charge  effects .  The  agreement  between  the  theory  for  space  charge 

limited  currents  and  the  data  presented  in  Fig.  18  is  evident  from  the  following 

38 

niunerical  calculations.  Lampert  relates  the  break  point  voltage  (that  is,  the 
voltage  at  which  space  charge  effects  begins  to  take  over  the  conduction  mechanism) 
V  (fi  -  s),  the  resistivity  p  ,  the  dielectric  constant  e  ,  the  mobility  M  and 
the  electrode  spacing  /  in  the  following  equation: 

g 

With  e  =  10,  p  =  1,  i  -  0.2  cm,  p=  4.4x10  fi-cm  (as  determined  from 

the  initial  slope  of  the  curve),  this  equation  yields  a  value  of  V(f2  -  s)  of  50 

volts.  The  experimental  value  is  44  volts .  Furthermore,  a  calculation  of  the 

39 

current  at  a  high  voltage,  according  to  an  equation  given  by  Bube  gives,  for 
the  total  current  J  , 
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-1*^  2 

10  V  eM  A  ^ 

J  =  - 5 — *- —  =  2  X  10  amp 

r 

for  a  voltage  of  200  volts  and  a  cross-sectional  area  A  of  4x  10”^  cm^.  The 

_9 

experimental  value  is  6x10  amp.  The  agreement  is  quite  good.  It  might 
appear  that  the  resistivity  of  the  sample  is  too  low.  Actually  it  was  necessary 
to  illuminate  the  sample  in  order  to  get  current  readings  within  the  limits  set 
by  the  electrometer  leakage  currents.  The  incident  radiation  was  345  m/x. 

Sample  J-5  might  have  shown  space  charge  limited  currents  if  a  high  enough 
field  had  been  applied . 

2.  Current  Saturation 

All  other  samples  (illuminated  or  in  dark  conditions)  displayed  current-voltage 
characteristics  of  the  type  shown  in  Fig.  19.  This  is  the  type  of  I-V  charac¬ 
teristic  one  would  expect  for  a  sample  with  the  type  of  band  structure  envisioned 
to  account  for  the  anomalous  photovoltaic  effect,  i.e. ,  opposing  barriers  at 
each  interface.  Thus,  the  crystal  should  act  like  a  number  of  pairs  of  diodes 
placed  back  to  back.  If  one  considers  just  two  identical  diodes  back-to-back, 
the  relation  between  the  current  (!)  and  the  voltage  (V)  across  both  diodes  is 
given  by  the  expression, 


1=1^  tanh  a  V  ,  (33) 

where  I  is  the  value  of  the  current  at  saturation,  a  is  e/2kT  for  two  identi- 
o 

cal  diodes.  The  results  of  plotting  I/I^  versus  tanh  /3V  is  presented  in  Fig. 

20.  I  was  estimated  to  be  0.050x10  amp  and  /3  was  chosen  as  0.0125. 
o  - 

The  points  generally  fall  on  straight  lines  but  the  slopes  of  the  two  lines  are  dif¬ 
ferent.  Here  again,  it  must  be  emphasized  that  the  junctions  in  the  crystals  are 
not  all  alike,  as  assumed  above. 
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Fig.  19  Current-Voltage  Curvee  for  a  Typical  Crystal  Showing  Saturation  Effect 
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*  Current-voltage  curves  for  a  typical  crystal  showing  saturation  effect.  The  two 
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Since  the  expression  for  the  current  is  symmetric  with  respect  to  voltage,  it  is 
clear  that  one  can  predict  the  behavior  of  any  number  of  diode  pairs  by  simply  re¬ 
placing  V  in  Equation  (33)  by  V/N,  where  N  is  just  the  number  of  diode  pairs. 
Conversely,  one  should  be  able  to  find  the  number  of  diode  pairs  from  the  experi¬ 
mental  data.  The  constant  /3,  chosen  above  to  give  a  good  fit  between  Equation 
(33)  and  the  data,  now  becomes  e/2NkT  where  N  is  the  number  of  diode  pairs. 
Using  this  argument,  it  was  found  that  these  crystals  generally  contain  about 
2,000  diode  pairs  per  cm.  According  to  the  model  of  the  APE,  this  means  that 
there  are  about  4,  000  layers/cm  in  the  crystal.  Merz,  ^  using  optical  techniques, 
arrived  at  a  figure  of  about  2 , 000  layers/cm . 
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VI.  CONCLUSIONS 

It  was  shown  that  earlier  models  failed  to  account  for  all  the  features  of  the 
APE.  Specifically,  both  Tanc's  model  and  Neumark's  model  could  not  account 
for  the  anomalously  high  voltages  at  both  the  positive  and  negative  peaks,  and 
the  effect  of  polarized  light  on  the  spectral  response  curve.  In  addition,  Tanc's 
scheme  to  add  elemental  voltages  proved  to  be  inadequate.  The  only  objection 
to  Hutson's  model  was  that  it  was  incomplete. 

By  extending  Hutson's  model  and  formulating  the  problem  in  terms  of  independent 
photovoltage  generators  in  each  segment  of  the  crystal,  one  can  account  for  all 
the  features  of  the  spectral  response  of  the  APE.  Briefly,  the  present  model 
correctly  predicts  the  shape  of  the  spectral  response  curve,  the  wavelengths  at 
which  the  photovoltage  changes  polarity,  and  the  general  position  of  the  photo¬ 
voltage  maxima.  Furthermore,  the  model  predicts  the  correct  dependence  of 
photocurrent  on  the  polarization  of  incident  light  and  the  shift  of  and 
with  temperature  and  light  intensity. 

The  expression  derived  for  the  dependence  of  the  photovoltage  on  an  external 
applied  stress,  together  with  the  results  of  the  stress  experiments  on  mixed 
crystals,  verifies  the  assumption  of  the  origin  of  the  photovoltage,  viz. ,  that 
the  photovoltages  arise  from  piezoelectric  fields  set  up  by  permanent  internal 
strains  formed  during  the  growth  process . 

The  experiment  on  single  crystal  hexagonal  ZnS  not  only  added  strength  to  the 
proposed  origin  of  the  APE ,  but  also  successfully  tied  together  the  information 
on  direction  of  crystal  growth,  pyroelectric  effect,  and  the  positive  peak  polarity. 
This  experiment  also  verified  the  calculation  of  the  photovoltage  generated  in  a 
non-uniformly  strained  piezoelectric  crystal. 
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In  general  then,  the  proposed  model  is  quite  successful;  the  main  weakness 
of  the  model  being  the  fact  that  the  photovoltage  depnds  on  the  internal  strain 
gradients ,  quantities  which  cannot  be  independently  measured. 
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THEVENIN'S  THEOREM  FOR  A  PHOTOVOLTAIC  GENERATOR 


22 


According  to  Tauc,  the  field  F  is  related  to  the  current  I,  the  conductivity 
cr  ,  tj  ,  and  t^  the  transfer  coefficients  for  electrons  and  holes,  and  the 
chemical  potentials  for  electrons  and  holes  £  ^  and  £  2 


F  =  i_il£i  +  !2lf2 
~  e  dx  e  dx 


To  find  the  short-circuit  current,  it  is  necessary  to  solve  the  equation  with  the 
condition  that 


/ 


Fdx  =  0 


Thus, 


At  =  i  f 

j  e^/yldx’Sdx^ 


dx  . 


But  since  I  is  constant  with  respect  to  x  ,  this  equation  can  be  rewritten  as. 


if(t  ^ 

e  *  \  1  dx 


d£. 

^2  1^) 


sc 


dx 


or  replacing  a  by  p  ,  Ig^,  becomes 


'sc  = 


if(t  ^ 

ejyi  dx 


d£, 


'2 
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22 

Again  referring  to  Tauc,  one  finds  that 


'  “OC  ■ 


where  Uq^  Is  the  open-circuit  voltage.  Thus, 


^8C  ~  ^f>  dx 

If  it  is  assumed  that  the  photogenerator  extends  from  x  =  a  to  x  =  b  ,  has 
a  unit  cross-section,  and  that  the  circuit  is  completed  by  a  resistanceless  wire, 
then 


sc 


PC 
p  dx 


But, 


where  R  is  the  total  resls^'>r,r’f>  between  a  and  b  in  the  light.  Thus,  I,.,  can  be 
written 


'sc  = 


= 

R 
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